METHOD AND APPARATUS FOR ESTIMATING POSITION UTILIZING GPS 

SATELLITE SIGNAL 



BACKGROUND OF THE INVENTION 

5 

This application claims the priority of Korean Patent Application No. 2003-0004579, 
filed on January 23, 2003, in the Korean Intellectual Property Office, the contents of 
which are incorporated herein in their entirety by reference. 

10 1. Field of the Invention 

The present invention relates to a method and apparatus by which the position 
of a position estimation apparatus is estimated using a plurality of satellite signals that 
are output with predetermined time differences fi:om one or more global positioning 
system (GPS) satellites. 

15 

2. Description of the Related Art 

Generally, position estimation methods using the GPS satellite system estimate 
a range between the GPS satellite and the antenna of a position estimation apparatus (for 
example, a GPS receiver), using triangulation and clean and acquisition (C/A, or coarse 

20 and acquisition) code which is output from the GPS satellite. 

The GPS satellite always transmits an LI frequency of 1575.42MHz that 
carries a C/A code. A position estimation apparatus generates the same code as the 
C/A code. The generated C/A code is compared with the received C/A code of the 
GPS satellite. From the result of the comparison, the time it takes for the C/A code 

25 output from the GPS satellite to arrive at the position estimation apparatus is measured. 
The position estimation apparatus measures the range between the GPS 
satellite and the position estimation apparatus by using the result of multiplication of the 
speed of light (the speed of the C/A code output firom the GPS satellite) and the time 
taken for arrival. Since the C/A code comprises a pseudo random noise code that is 

30 almost noise in itself, and the measured range between the GPS satellite and the position 
estimation apparatus contains errors, the range is referred to as a pseudo range. 
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A conventional position estimation apparatus receives satellite signals 
simultaneously output from at least four GPS satellites, measures the pseudo range 
corresponding to each satellite signal, and from the measured results estimates the 
position of the position estimation apparatus. That is, the position estimation apparatus 
5 estimates a 3-dimensional position of the position estimation apparatus by 

simultaneously using four or more satellite signals output from four or more GPS 
satellites. 

However, due to the influence of the environment, the position estimation 
apparatus at times cannot receive at least four signals at the same time. In this case, 
10 the conventional position estimation apparatus cannot estimate its 3-dimensional 
position. 

SUMMARY OF THE INVENTION 

The invention is directed to a method and apparatus for determining position 

15 using a global position satellite (GPS) signal. A firist GPS signal is received at a 
receiver from a first GPS satellite having a first position. A second GPS signal is 
received from the satellite at a second position of the satellite. The position of the 
receiver is determined using the first and second GPS signals. 

A third GPS signal can be received from the first GPS satellite at a third 

20 position of the first GPS satellite, and the third GPS signal can be used in determining 
the position'of the receiver. A fourth GPS signal can be received from the first GPS 
satellite at a fourth position of the first GPS satellite, and the fourth GPS signal can be 
used in determining position of the receiver. The fourth GPS signal can also be 
received from a second GPS satellite at a first position of the second GPS satellite, and 

25 this fourth GPS signal can be used in determining the position of the receiver. 

The third GPS signal can also be received from a second GPS satellite at a first 
position of the second GPS satellite. This third GPS signal can be used in determining 
the position of the receiver. The fourth GPS signal can be received from the second 
GPS satellite at a second position of the second GPS satellite, and this fourth GPS signal 

30 can be used in determining position of the receiver. This fourth GPS signal can also be 
received from a third GPS satellite at a first position of the GPS satelHte, and this fourth 
GPS signal can be used in determining position of the receiver. 
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The apparatus for determining position includes a receiver for receiving the 
GPS signals and a position calculation unit for determining position of the receiver. A 
controller detects a number of usable satellites that can be used in determining position. 
A stationary measurement request and selection unit can request that the user remain 
5 stationary during determination of position. The position calculation unit can include a 
time difference measurement determiner which requests that the user remain stationary 
during determination of position if a number of usable satellites is below a threshold. 
The position calculation unit also includes a time difference measurement calculator 
which calculates position by measuring time differences between the GPS signals, 
10 In accordance with the invention, position of a receiver such as a GPS receiver 

can be determined using fewer than four GPS signals from GPS satellites. As a result, 
a much more reliable position determination or estimation system is provided, 
especially in cases where conditions such as environmental conditions prevent reception 
of GPS signals from four or more satellites. 

15 

BRIEF DESCRIPTION OF THE DRAWINGS 

The foregoing and other objects, features and advantages of the invention will 
be apparent from the following more particular description of preferred embodiments of 
20 the invention, as illustrated in the accompanying drawings in which hke reference 
characters refer to the same parts throughout the different views. The drawings are not 
necessarily to scale, emphasis instead being placed upon illustrating the principles of the 
invention. 

FIG 1 is a diagram illustrating a position estimation method according to a first 
25 embodiment of the present invention. 

FIG 2 is a diagram illustrating in more detail the position estimation method 
according to the first embodiment of the present invention. 

FIG. 3 is a diagram illustrating a position estimation method according to a 
second embodiment of the present invention. 
30 FIG 4 is a diagram illustrating in more detail the position estimation method 

according to the first embodiment of the present invention. 

FIG S is a schematic block diagram of a position estimation apparatus 
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according to an embodiment of the present invention. 

FIG 6 is a detailed schematic block diagram of the controller and the position 
calculation unit shown in FIG 5. 

FIG 7 is a flowchart illustrating steps performed by a position estimation 
5 method based on the number of satellites that can be measured according to the present 
invention. 

FIG. 8 is a flowchart illustrating in more detail the position estimation method 
of a position estimation apparatus as shown in FIG 5 by measuring time differences, in 
accordance with the present invention. 

10 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Referring to FIG 1, reference number 11 indicates the position of a GPS 
satellite at time (tl), and reference number 12 indicates the position of the GPS satellite 
at time (t2) after a lapse of a predetermined period from time (tl). 

1 5 Accordingly, if each of the pseudo ranges of the satellite signals which are 

output at time (tl) and time (t2) from the identical GPS satellite is measured and then 
results of the two range equations are obtained, the estimated positions of the position 
estimation apparatus, for example a GPS receiver, may be two places (13 A, 13B) if the 
positions are in a plane. 

20 However, since the position of the position estimation apparatus is in space 

coordinates and is variable, and there is a time difference between the GPS satellite 
measurements for the position estimation apparatus, four range equations are required in 
order to actually obtain the position of the position estimation apparatus. But, the 
method for positioning the position estimation apparatus is not limited to the method 

25 using the pseudo range. 

FIG 2 is a diagram illustrating in more detail the position estimation method 
according to the first preferred embodiment of ttie present invention. In the position 
estimation method according to the present invention, a plurality of satellite output 
signals with predetermined time differences from one or more GPS satellites are 

30 received, and in response to each of the received satellite signals, the position of the 
position estimation apparatus is estimated. 

Referring to FIG.2, coordinates (x, y, z) indicate an estimated position of a 



SAM-0444 



position estimation apparatus 25, coordinates (xl, yl, zl, tl) indicate position data (xl, 
yl, zl) of a GPS satellite 21 at time (tl), and pi indicates the pseudo range measured at 
time (tl) by the position estimation apparatus. Here, it is preferable that the position 
estimation apparatus 25 be stationary or fixed. 

5 Coordinates (x2, y2, z2, t2) indicate position data (x2, y2, z2) of a GPS satellite 

22 at time (t2), and p2 indicates the pseudo range measured at time (t2) by the position 
estimation apparatus. Here, time (t2) indicates a lapse of a predetermined period from 
time (tl). 

Coordinates (x3, y3, z3, t3) indicate position data (x3, y3, z3) of a GPS satellite 
10 23 at time (t3), and p3 indicates the pseudo range measured at time (t3) by the position 
estimation apparatus. Here, time (t3) indicates a lapse of a predetermined period from 
time (t2). 

Coordinates (x4, y4, z4, t4) indicate position data (x4, y4, z4) of a GPS satellite 
24 at time (t4), and p4 indicates the pseudo range measured at time (t4) by the position 
15 estimation apparatus. Here, time (t4) indicates a lapse of a predetermined period from 
time (t3). Reference numbers 21-24 do not indicate four different satellites, but the 
same satellite at four different times and locations. 

The identification number (ID) of each of the GPS satellites 21 through 24 may 
be identical or different numbers. For example, the IDs of at least two or more GPS 
20 satellites may be identical numbers. 

The following equation 1 shows four range equations expressed by respective 
position data ((xl, yl, zl), (x2, y2, z2), (x3, y3, z3), (x4, y4, z4)) of the GPS satellites, 
and measured pseudo ranges (pi, p2, p3, p4) at points of time (tl, t2, t3, t4). 

Accordingly, the position estimation apparatus receives position data ((xl, yl, 
25 zl), (x2, y2, z2), (x3, y3, z3), (x4, y4, z4)) that are output from respective GPS satellites 
21 through 24 at respective points of time (tl, t2, t3, t4). Then, by measuring 
respective pseudo ranges (pi, p2, p3, p4) and obtaining solutions (or roots) for the four 
range equations given as equation 1, the position 25, for example, space coordinates (x, 
y, z), of the position estimation apparatus can be calculated or estimated. 

30 
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^{xl - xf + (yl - yf + (zl - + cAt = pi 
yl{x2 - xY + (y2 - yf + (z2 - zf + cAt = p2 
^|(x3 - xf + iy2 - yf + (z3 - zf ^cAt = p3 
V(a:4 - jc)' + (y4 - yf + (z4 - z)' + cA/ = p4 (1) 

Here, c denotes the speed of light and At denotes the difiference between the 
time of a GPS satellite and the time of the position estimation apparatus. 
5 In order to obtain space coordinates (x, y, z) 25 of the position estimation 

apparatus, at least four range equations as shown in equation 1 are needed. 
Accordingly, the number of GPS satellites whose satellite signals can be simultaneously 
received determines the minimum number of measuring times of pseudo ranges. 

For example, if the position estimation apparatus can receive satellite signals 
10 output from two GPS satellites, the position estimation apparatus receives the satellite 
signals output from the two GPS satellites at a point of time and measures respective 
pseudo ranges. Then, at another point of time, the position estimation apparatus 
receives the satellite signals output from the two GPS satellites (or two satellite signals 
output from two GPS satellites) and measures respective pseudo ranges. By doing so, 
15 four range equations as shown in equation 1 can be obtained. 

In this case, in order to obtain the position 25 of the position estimation 
apparatus, that is, the space coordinates (x, y, z), the pseudo range is measured at least 
twice at each of the different points of time. 

Also, if the position estimation apparatus can receive satellite signals output 
20 from three GPS satellites, the position estimation apparatus receives the satellite signals 
output from the three GPS satellites at a point of time and measures respective pseudo 
ranges. Then, at another point of time, the position estimation apparatus receives the 
satellite signal output from any one of the thiree GPS satellites and measures the pseudo 
range. By doing so, four range equations as shown in equation 1 can be obtained. 
25 Accordingly, in order to obtain the position 25 of the position estimation 

apparatus, that is, the space coordinates (x, y, z), the pseudo range is measured at least 
twice at each of the different points of time. 

If the altitude (that is, the coordinate of the z-axis) of the position estimation 

SAM-0444 R 



apparatus is already known, the minimum number of measuring times of pseudo ranges 
for estimating the space coordinates of the position estimation apparatus is less by at 
least once than that for estimating the space coordinates when the altitude is not known. 
FIG 3 is a conceptual diagram illustrating a position estimation method 
5 according to a second preferred embodiment of the present invention. The position 
estimation apparatus estimates the position of the position estimation apparatus, by 
using the differences of every two pseudo ranges that are close in time. 

In this method, a pseudo range is measured from each of a plurality of satellite 
signals that are output at relatively shorter time differences from a single GPS satellite 

10 (for example, the ID of the GPS satellite, each point of time, position data at respective 
points of time). Then, by using the measured pseudo ranges, the position of the 
position estimation apparatus is estimated. Accordingly, delays common to the 
satellite signals, for example, ionospheric and tropospheric delay, cancel each other such 
that the position of the position estimation apparatus can be more precisely estimated. 

15 When the position of the position estimation apparatus is estimated by using 

the differences of every two pseudo ranges that are close in time, the estimated positions 
of the position estimation apparatus correspond to coordinates on hyperbolas. 
Accordingly, intersections where two or more hyperbolas intersect are the estimated 
positions of the position estimation apparatus. 

20 Referring to FIG 3, reference numbers 3 1, 32, and 33 indicate positions of a 

single GPS satellite, which has an identical ID at different points of time, and pi, p2, 
and p3 indicate pseudo ranges measured by the position estimation apparatus at 
positions 31, 32, and 33, respectively. 

In a plane, hyperbola (kl), in which the value of (pl-p2) is constant, and 

25 hyperbola (k2), in which the value of (p2-p3) is constant, meet on two points (34A, 
34B) such that the two points (34A, 34B) correspond to an estimated position of the 
position estimation apparatus. 

However, since flie position of the position estimation apparatus is a 3- 
dimensional space position, the pseudo ranges must be measured at least four times at 

30 different points of time in order to estimate the position of the position estimation 
apparatus. If the altitude (that is, the coordinate of the z-axis) of the position 
estimation apparatus is known, the position estimation apparatus can estimate the spatial 
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position of the position estimation apparatus by measuring pseudo ranges at least three 
times at different pomts of time. 

FIG 4 is a diagram illustrating in more detail the position estimation method 
according to the first preferred embodiment of the present invention. Coordinates (x, y, 
5 z) indicate an estimated position of a position estimation apparatus 45, coordinates (xl, 
yl, zl, tl) indicate position data (xl, yl, zl) of a GPS satellite 41 at time (tl), and pi 
indicates the pseudo range measured at time (tl) by the position estimation apparatus. 
Coordinates (x2, y2, z2, t2) indicate position data (x2, y2, z2) of a GPS satellite 42 at 
time (t2), and p2 indicates the pseudo range measured at time (t2) by the position 

10 estimation apparatus. Here, time (t2) indicates a lapse of a predetermined period from 
time (tl). Coordinates (x3, y3, z3, t3) indicate position data (x3, y3, z3) of a GPS 
satellite 43 at time (t3), and p3 indicates the pseudo range measured at time (t3) by the 
position estimation apparatus. Here, time (t3) indicates a lapse of a predetermined 
period from time (t2). Coordinates (x4, y4, z4, t4) indicate position data (x4, y4, z4) of a 

15 GPS satellite 44 at time (t4), and p4 indicates the pseudo range measured at time (t4) by 
the position estimation apparatus. Here, time (t4) indicates a lapse of a predetermined 
period from time (t3). Here, the GPS satellites 41 through 44 are the same GPS 
satellite, and it is preferable that the position estimation apparatus 45 be stationary or 
fixed. 

20 The following equation 2 shows four pseudo range equations expressed by 

respective position data ((xl, yl, zl), (x2, y2, z2), (x3, y3, z3), (x4, y4, z4)) of the GPS 
satellite at points of time (tl, t2, t3, t4), and measured pseudo ranges (pi, p2, p3, p4). 
Here, kl denotes a hyperbola in which the value of (pl-p2) is constant, k2 denotes a 
hyperbola in which the value of (p2-p3) is constant, and k3 denotes a hyperbola in 

25 which the value of (p3-p4) is constant. 

Accordingly, the position estimation apparatus receives position data ((xl, yl, 
zl), (x2, y2, z2), (x3, y3, z3), (x4, y4, z4)) that are output from the same GPS satellite at 
respective points of time (tl, t2, t3, t4). Then, by measuring respective pseudo ranges 
(pi, p2, p3, p4) and obtaining solutions (or roots) for the range equations given as 

30 equation 2, position 45, for example, space coordinates (x, y, z), of the position 
estimation apparatus can be calculated or estimated. 
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-xY+iyl-yy+izl-zf- 4{x2 - xf ^iyl-yf^ (z2 - zf = k\ 



4(x2 -xf^iyl-yY-^ (z2 - zf - ^(^3 -xY^{y3- yf + (z3 - zf = kl 
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FIG 5 is a schematic block diagram of a position estimation apparatus 
according to a preferred embodiment of the present invention. Referring to FIG 5, the 
5 position estimation apparatus (for example, a GPS receiver 500) comprises an antenna 
510, a signal processing unit 520, a position calculation unit 540, and a stationary 
measurement request and selection unit 550. 

The antenna 510 receives a satellite signal from a satellite. The signal 
processing unit 520 calculates correlation values between C/A codes loaded on the 
10 satellite signal received through the antenna 510 and self-generated C/A codes and then 
outputs delay information based on the calculation result to the position calculation unit 
540. 

The signal processing unit 520 includes a pre-amplifier 521, a down converter 
523, an analog-to-digital (A/D) converter 525, an automatic gain controller (AGC) 527, 
1 5 a mixer 529, a carrier numerically controlled oscillator (NCO) 53 1 , a code generator 
533, a correlator 535, a code NCO 537 and a controller 539. 

The pre-amplifier 521 amplifies the satellite signal received from the satellite 
through the antenna 510 and outputs the amplified satellite signal to the down converter 
523. 

20 The down converter 523 receives the amplified satellite signal, converts the 

same into an intermediate frequency signal in response to the output signal of the AGC 
527, and outputs the converted signal to the A/D converter 525. 

The A/D converter 525 receives the intermediate frequency signal and converts 
the same into a digital signal. 
25 The AGC 527 controls the gain of the down converter 523 in response to the 

output signal of the A/D converter 525. 

The carrier NCO 531 generates I (in phase) sine waves and a Q (quadrature 
phase) sine waves in order to compensate for a doppler effect of the satellite signal. 

The mixer 529 mixes the output signal of the A/D converter 525 with the I and 
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Q sine waves output from the carrier NCO 531 to offset the doppler effect of the 
satellite signal. Thus, the mixer 529 provides C/A codes output from a GPS satellite to 
the correlator 535. 

The code NCO 537 generates delayed codes in accordance with an expected 
5 delay of the satellite signal from the C/A codes generated from the code generator 533, 
and outputs the generated codes to the correlator 535. The code generator 533 
generates C/A codes based on a reference time of the position estimation apparatus 500 
and an identification (ID) number of a satellite to be detected. 

The correlator 535 calculates a correlation value between the C/A code loaded 
10 on the output signal of the mixer 529 and the C/A code output from the code NCO 537, 
and outputs the delay information based on the calculation result to the controller 539. 

The controller 539 detects the satellite sfignal in response to the output signal of 
the correlator 535. If no satellite signal is detected, the controller 539 controls the next 
expected doppler frequency, delayed code value or a code value of another satellite to 
1 5 be output to the mixer 529 or the correlator 535 through the carrier NCO 53 1 and the 
code NCO 537. Thus, the operation of the correlator 535 is iteratively performed. 

However, if a desired satellite signal is detected, the controller 539 controls the 
calculation result of the correlator 535 to be output to the position calculation unit 540. 
The position calculation unit 540, in response to the delay signal output from 
20 the signal processing unit 520, obtains pseudo ranges and/or solutions of distance 
equations expressed by equations 1 and 2. In this case, when the number of 
simultaneously measurable satellites is generally smaller than that of satellites required 
for estimating the position of the position estimation apparatus 500, the position 
calculation unit 540 requests a user of the position estimation apparatus 500 to stop or 
25 remain stationary by voice or text through the stationary measurement request and 
selection unit 550 using the position estimation method according to the present 
invention. 

After the user stops, the stationary measurement request and selection unit 550 
selects a stationary measuring function to then output a predetermined selection signal 
30 to the position calculation unit 540. Thus, the position calculation unit 540 starts 

estimation of the stationary position of the position estimation apparatus 500 using the 
position estimation method according to the present invention. 
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FIG 6 is a detailed diagram of the controller 539, the stationary measurement 
request and selection unit 550 and the position calculation unit 540 shown in FIG 5. 
Referring to FIGS. 5 and 6, the stationary measurement request and selection unit 550 
includes a stop request display 5501 and a stop selector 5503. The controller 539 
5 includes a stationary signal processor 5391 and a receiver controller 5393. The 

position calculation unit 540 includes a time difference measurement determiner 5401, a 
time difference measurement calculator 5403 and a GPS position calculator 5405. 

The time difference measurement determiner 5401 determines whether time 
difference measurement is initiated or not, in response to the delay information output 

10 from the controller 539. The time difference measurement calculator 5403 calculates 
the position of the position estimation apparatus 500 by measuring the time difference. 
The GPS position calculator 5405 calculates the position of the position estimation 
apparatus 500 by an ordinary method. That is, the GPS position calculator 5405 is an 
exemplary circuit for estimating the position of the position estimation apparatus 500 

15 from a simultaneously received plurality of satellite signals. 

If the receiver controller 5393 outputs delay information corresponding to data 
of detected satellites, that is, satellite signals output from measurable satellites, to the 
time difference measurement determiner 5401 of the position calculation unit 540, the 
time difference measurement determiner 5401 determines whether there enough satellite 

20 signals to calculate the position of the position estimation apparatus 500, based on the 
received delay information. 

If the number of received satellite signals is sufficient, the time difference 
measurement determiner 5401 instructs the general GPS position calculator 5405 to 
calculate the position of the position estimation apparatus 500. If the number of 

25 received satellite signals is not sufficient, the time difference measurement determiner 
5401 sends a request signal for stationary measurement to the stationary signal 
processor 5391. 

The stationary signal processor 5391 which has received the signal output from 
the time difference measurement determiner 5401 notifies the user to stop through the 
30 stop request display 5501 . Then, if the notified user permits stationary measurement 
through the stop selector 5503, the stop selector 5503 outputs a permission signal to the 
time difference measurement determiner 5401 through the stationary signal processor 
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5391. 

The time difference measurement determiner 5401 having received the 
permission signal instructs the time difference measurement calculator 5403 to start 
time difference measurement according to the present invention on the assumption that 
5 the user has stopped. 

When the stationary m'easurement is not permitted even if the user has been 
requested to stop through the stop request display 5501, the time difference 
measurement determiner 5401 cannot calculate the position of the position estimation 
apparatus 500 and continuously detects satellites necessary for calculating the position 
1 0 of the position estimation apparatus 500. 

Since the user remains stationary in the course of estimating the position of ttie 
position estimation apparatus 500 using a time difference, the user is continuously 
notified to remain stationary. For example, in order to start time difference measurement, 
the stop request display 5501 flickers a light. The light is on while measuring a time 
1 5 difference. Thus, the user can distinguish a stop request for starting time difference 
measurement and a stop request for maintaining the stationary state. 

Also, if time difference measurement is not necessary due to reception of 
enough satellite signals while measuring a time difference, the stop request display 5501 
turns the light off, thereby notifying the user that he/she need not remain stationary. 
20 FIG. 7 s a flowchart of steps performed in a position estimation method based 

on the number of satellites that can be measured, of the first and second preferred 
embodiments of the present invention. Referring to FIGS. 5 through 7, the position 
calculation unit 540, in response to a plurality of satelUte signals that are output with 
predetermined time differences from at least one or more GPS satellites, can estimate 
25 the position of the position estimation apparatus 500. 

The position calculation unit 540 receives a plurality of satellite signals that are 
output with predetermined time differences from one or more GPS satellites and then 
can estimate the position of the position estimation apparatus 500. 

If the position estimation apparatus 500 begins position estimation in step 70, 
30 then the time difference measurement determiner 5401 of the position calculation unit 
540 determines whether or not the number of satellites whose pseudo range can be 
measured (hereinafter referred to as "the number of measurable satellites") is greater 
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than 3 in step 71. That is, if satelUte signals from four or more satellites can all be 
received, or if the number of measurable satellites is greater than 3, the general GPS 
position calculation unit S40S of the position calculation unit 540 estimates its position 
in the ordinary method in step 72. 
5 When the number of measurable satellites is 3 in step 73, if the altitude 

information (the coordinate z) of the position estimation apparatus 500 is already known 
in step 74, the general GPS position calculation unit 5405 of the position calculation 
unit 540 estimates its position in the ordinary method in step 75. 

However, when the number of measurable satellites is 1 or 2 in step 76, 

10 or when the number of measurable satellites is 3 and the altitude information of the 
position estimation apparatus 500 is not known in steps 73 and 74, the stop request 
display 5501 asks the user of the position estimation apparatus 500 to stop and estimates 
the position of the position estimation apparatus 500 according the position estimation 
method of the present invention. If the user, in response to this, stops and begins to 

15 execute the position estimation method according to the present invention in step 77, the 
position estimation apparatus 500 measures a pseudo range in response to each of a 
plurality of satellite signals that are input with predetermined time differences 
(hereinafter referred to as "time difference measurement"), and estimates the position of 
the position estimation apparatus 500 based on the satellite signals and the respective 

20 pseudo ranges in step 78. In this case, equation 2 is used. 

However, if there is no measurable satellites at all in step 76, or if the user does 
not respond to the stop request and does not execute the position estimation method in 
step 77, it is not possible to estimate the position of the position estimation apparatus 
500 in step 79. 

25 FIG 8 is a flowchart showing in more detail the position estimation method 

according to the first and second embodiments of the present invention. That is, FIG 8 
is a flowchart illustrating step 78 of FIG 7 in more detail. 

If the position estimation of the position estimation apparatus 500 by the time 
difference measurement begins in step 78, the time difference measurement determiner 

30 5401 of the position calculation unit 540 determines whether or not the number of 
measurable satellites is 3 in step 81. If the determination result indicates that the 
number of measurable satellites is 3 in step 81, the position estimation apparatus 500 
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measures a pseudo range in response to each of satellite signals that are output from the 
three measurable satellites at a first point of time, and measures a pseudo range in 
response to a satellite signal that is output from any one of the three measurable 
satellites at a second point of time after a lapse of predetermined time from the first 

5 point of time, or a satellite signal output from a new satellite other than the three 

satellites, in step 83. Accordingly, the minimum number of measuring time differences 
after the first point of time is one. 

If the number of measurable satellites is 2 in step 84, the time difference 
calculator 5403 measures a pseudo range in response to each of satellite signals that are 

10 output form the two measurable satellites at a first point of time, and measures a pseudo 
range in response to each of satellites signals tiiat are output from the two measurable 
satellites at a second point of time after a lapse of predetermined time from the first 
point of time, or a satellite signal output from a new satellite other than ttie two satellites, 
in step 85. Accordingly, the minimum number of measuring time differences after the 

1 5 first point of time is two. 

If the number of measurable satellites is 2 in step 84, and the altitude 
information of the position estimation apparatus 500 is already known, the position 
estimation apparatus 500 measures a pseudo range in response to each of the satellite 
signals output from the two measurable satellites at the second point of time in step 85. 

20 Accordingly, the minimum number of measuring time differences after the first point of 
time is one. 

If the number of measurable satellites is not 2, that is, if the number of 
measurable satellites is 1 in step 84, the position estimation apparatus 500 measures a 
pseudo range in response to the satellite signal output from the measurable satellite at a 

25 first point of time, measures a pseudo range in response to the satellite signal output 
from the measurable satellite at a second point of time after a lapse of predetermined 
time from the first point of time, or a satellite signal output from a new satellite other 
than the satellite, measures a pseudo range in response to the satellite signal output from 
the measurable satellite at a third point of time after a lapse of predetermined time from 

30 the second point of time, or a satellite signal output from a new measurable satellite, and 
measures a pseudo range in response to the satellite signal output from the measurable 
satellite at a fourth point of time after a lapse of predetermined time from the third point 
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of time, or a satellite signal output from a new measurable satellite. Accordingly, the 
minimum number of measuring time differences after the first point of time is 3 times in 
step 86. 

However, if the altitude information of the position estimation apparatus 500 is 
5 already known, the position estimation apparatus 500 measures a pseudo range in 
response to the satellite signal output from the measurable satellite at the second point 
of time, or a satellite signal output from a new measurable satellite, and measures a 
pseudo range in response to the satellite signal output from the measurable satellite at 
the third point of time, or a satellite signal output from a new measurable satellite. 
1 0 Accordingly, the minimum number of measuring time differences after the first point of 
time is two in step 86. 

When the time difference measurement is performed, the position estimation 
apparatus 500 determines at each point of time whether or not a measurable satellites is 
the same satellite in step 87. If the determination result indicates that the satellites is 
15 the same satellite, the position estimation apparatus 500 measures the pseudo range 

between the position estimation apparatus 500 and the satellite with the time difference, 
and estimates the position of the position estimation apparatus 500 by using the 
differences of every two pseudo ranges expressed by equation 2 in step 88. 

If the determination result indicates that the satellites are not an identical 
20 satellite, the position estimation apparatus 500 estimates its position by measuring the 
pseudo range between the position estimation apparatus 500 and the satellite with the 
time difference in step 89. 

As described above, by the position estimation method and apparatus 
according to the present invention, even when the number of measurable satellites is 3 
25 or less, the pseudo ranges can be measured by using time differences, and by using the 
measured pseudo ranges, the position of the position estimation apparatus can be 
precisely estimated or calculated. 

While this invention has been particularly shown and described with reference 
to preferred embodiments thereof, it will be understood by those skilled in the art that 
30 various changes in form and details may be made therein without departing from the 
spirit and scope of the invention as defined by the appended claims. 
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